source("~/jules_ppe_gmd/docs/JULES-ES-1p0-common-packages.R")
## Loading required package: spam
## Loading required package: dotCall64
## Loading required package: grid
## Spam version 2.7-0 (2021-06-25) is loaded.
## Type 'help( Spam)' or 'demo( spam)' for a short introduction
## and overview of this package.
## Help for individual functions is also obtained by adding the
## suffix '.spam' to the function name, e.g. 'help( chol.spam)'.
##
## Attaching package: 'spam'
## The following objects are masked from 'package:base':
##
## backsolve, forwardsolve
## Loading required package: viridis
## Loading required package: viridisLite
## See https://github.com/NCAR/Fields for
## an extensive vignette, other supplements and source code
## Loading required package: lhs
## Loading required package: parallel
## Loading required package: viztools
source("~/jules_ppe_gmd/docs/JULES-ES-1p0-common-functions.R")
## ----------------------------------------------------------------------
## Data locations and constants
## ----------------------------------------------------------------------
#ensloc <- '/project/carbon_ppe/JULES-ES-1p0_PPE/'
ensloc_wave00 <- '/data/users/hadaw/JULES_ES_PPE/u-au932/'
ensloc_wave01 <- '/data/users/hadaw/JULES_ES_PPE/u-ck006/'
# Some pallete options
yg = brewer.pal(9, "YlGn")
ryb = brewer.pal(11, "RdYlBu")
byr = rev(ryb)
rb = brewer.pal(11, "RdBu")
br = rev(rb)
blues = brewer.pal(9, 'Blues')
cbPal <- c("#000000", "#E69F00", "#56B4E9", "#009E73", "#F0E442", "#0072B2", "#D55E00", "#CC79A7")
wave00col <- 'skyblue2'
wave01col <- 'tomato2'
zissou5 <- wes_palette('Zissou1', 5, type = c('discrete', 'continuous'))
zblue <- makeTransparent(as.character(zissou5)[1], 150)
zred <- makeTransparent(as.character(zissou5)[5], 150)
ysec = 60*60*24*365
years <- 1850:2013
global_mean_modern_value_wave00_file <- "~/jules_ppe_gmd/docs/data/global_mean_modern_value_wave00_2022-13-09.rdata"
load(global_mean_modern_value_wave00_file)
modern_value_stan_file <- "~/jules_ppe_gmd/docs/data/modern_value_stan_2022-09-13.rdata"
load(modern_value_stan_file)
Y <- datmat_wave00
Y_nlevel0_ix <- which(is.na(datmat_wave00[,'year']))
# Load up the data
lhs_i = read.table('~/jules_ppe_gmd/docs/data/lhs_u-ao732.txt', header = TRUE)
lhs_ii = read.table('~/jules_ppe_gmd/docs/data/lhs_u-ao732a.txt', header = TRUE)
toplevel_ix = 1:499
# The raw input data is a latin hypercube
lhs = rbind(lhs_i, lhs_ii)[toplevel_ix, ]
lhs_level0 <- lhs[-Y_nlevel0_ix,]
X = normalize(lhs)
colnames(X) = colnames(lhs)
X_level0 <- X[-Y_nlevel0_ix,]
X_nlevel0 <- X[Y_nlevel0_ix,]
d = ncol(X)
# lower and higher bound on the normalised matrix for visualisation
rx = rbind(rep(0,32), rep(1,32))
low_npp_ix <- which(Y[,'npp_nlim_lnd_sum'] < 1e5)
# code from https://stackoverflow.com/questions/28182872/how-to-use-different-sets-of-data-in-lower-and-upper-panel-of-pairs-function-in
#X <- matrix(runif(300), ncol=3)
#Y <- matrix(c(sort(runif(100, 0, 10)),
# sort(runif(100, 0, 10)),
# sort(runif(100, 0, 10))), ncol=3)
#pdf(file = 'figs/fig02.pdf', width = 12, height = 10)
#pdf(file = 'figs/run-failure-pairs.pdf', width = 12, height = 10)
x1 <- X[low_npp_ix, ]
x2 <- X_nlevel0
XY <- rbind(x1, x2)
pairs(XY,
lower.panel=function(x, y, ...) {
Xx <- x[seq_len(nrow(x1))] # corresponds to X subset
Xy <- y[seq_len(nrow(x1))] # corresponds to X subset
#usr <- par("usr"); on.exit(par(usr))
#par(usr = c(range(x1[, -ncol(x1)]), range(x1[, -1]))) # set up limits
points(Xx, Xy, col = zblue, pch = 19, cex = 0.8)
# if(par('mfg')[2] == 1) axis(2) # if left plot, add left axis
#if(par('mfg')[1] == ncol(x1)) axis(1) # if bottom plot add bottom axis
},
upper.panel=function(x, y, ...) {
Yx <- x[(nrow(x1) + 1):length(x)] # Y subset
Yy <- y[(nrow(x1) + 1):length(y)] # Y subset
#cntr <- outer(Yx, Yx, FUN='*') # arbitrary function for contour
# usr <- par("usr"); on.exit(par(usr))
#par(usr = c(range(x2[, -1]), range(x2[, -ncol(x2)]))) # set up limits
points(Yx, Yy, col = zred, pch = 19, cex = 0.8)
#contour(Yx, Yy, cntr, add=TRUE)
#if(par('mfg')[2] == ncol(x2)) axis(4) # if right plot, add right axis
#if(par('mfg')[1] == 1) axis(3) # if top plot, add top axis
},
#tick=FALSE, # suppress the default tick marks
#line=1,
gap = 0,
xlim = c(0,1), ylim = c(0,1),
labels = 1:d,
oma = c(2, 18, 2, 2)) # move the default tick labels off the plot
reset()
legend('left', legend = paste(1:d, colnames(lhs)), cex = 1.1, bty = 'n')
legend('topleft', pch = 19, col = c( zred, zblue), legend = c('failed', 'zero carbon cycle'), bty = 'n', inset = 0.02, cex = 1.1 )
#dev.off()
Y_char <- rep('RAN', nrow(X))
Y_char[Y_nlevel0_ix] <- 'CRASHED'
Y_char[low_npp_ix] <- 'LOWCARBON'
Y_factor <- as.factor(Y_char)
wave0_summary_df <- cbind(X, Y, Y_char, Y_factor)
save(lhs, X, Y, Y_char, Y_factor, wave0_summary_df, file = 'data/wave0_summary.RData')
Have we got the factors right?
cbind(Y[, 'npp_nlim_lnd_sum'], Y_char)
## Y_char
## [1,] "2464953.47849605" "RAN"
## [2,] "1248000.44246179" "RAN"
## [3,] "1893850.67745132" "RAN"
## [4,] "1779634.61879458" "RAN"
## [5,] "1333052.53686655" "RAN"
## [6,] "1829087.04066031" "RAN"
## [7,] "125.100277029329" "LOWCARBON"
## [8,] "1354509.40828833" "RAN"
## [9,] "1436905.38963356" "RAN"
## [10,] "1105586.84200311" "RAN"
## [11,] "2281524.08799275" "RAN"
## [12,] "1569806.15413172" "RAN"
## [13,] "1452198.55278638" "RAN"
## [14,] "1199097.02939196" "RAN"
## [15,] "2617350.39816265" "RAN"
## [16,] "1249447.10148417" "RAN"
## [17,] "939242.694289745" "RAN"
## [18,] "5.37100446251452" "LOWCARBON"
## [19,] "1077579.05514713" "RAN"
## [20,] "1523180.02583773" "RAN"
## [21,] "1569417.63070336" "RAN"
## [22,] "1873231.77735771" "RAN"
## [23,] "1912583.54670419" "RAN"
## [24,] "1902959.37750349" "RAN"
## [25,] "1445375.9153619" "RAN"
## [26,] "951627.465373486" "RAN"
## [27,] "1662356.03945454" "RAN"
## [28,] "1936958.91626488" "RAN"
## [29,] "1433603.10792285" "RAN"
## [30,] "976651.168260095" "RAN"
## [31,] "2082565.02331613" "RAN"
## [32,] "1531283.77037663" "RAN"
## [33,] "1645926.53544732" "RAN"
## [34,] "1320236.28995648" "RAN"
## [35,] "2891640.62239662" "RAN"
## [36,] "1635834.14721978" "RAN"
## [37,] "2089155.84530289" "RAN"
## [38,] "1377545.60029833" "RAN"
## [39,] "1487813.03878515" "RAN"
## [40,] "651768.709997753" "RAN"
## [41,] "2027624.59401214" "RAN"
## [42,] "2078839.96826463" "RAN"
## [43,] "518555.798220461" "RAN"
## [44,] "1061709.94504286" "RAN"
## [45,] "1775343.21596816" "RAN"
## [46,] "2311775.26816144" "RAN"
## [47,] "1441874.46802137" "RAN"
## [48,] "2195581.99198865" "RAN"
## [49,] "1325496.4911125" "RAN"
## [50,] "2004112.56580796" "RAN"
## [51,] "2368885.55711663" "RAN"
## [52,] "1361.54882660815" "LOWCARBON"
## [53,] "908731.084399189" "RAN"
## [54,] "1442169.69546151" "RAN"
## [55,] "139.629393179477" "LOWCARBON"
## [56,] "1333927.55011543" "RAN"
## [57,] "1666849.46116248" "RAN"
## [58,] "965867.767068223" "RAN"
## [59,] "1243866.85654177" "RAN"
## [60,] "1943266.73141755" "RAN"
## [61,] "1902138.05825955" "RAN"
## [62,] "1767409.94121306" "RAN"
## [63,] "1317749.66685722" "RAN"
## [64,] "1199841.24949906" "RAN"
## [65,] "1837154.16360185" "RAN"
## [66,] "1163145.1259366" "RAN"
## [67,] "1214317.39561894" "RAN"
## [68,] "1313614.75168373" "RAN"
## [69,] "1724789.67747131" "RAN"
## [70,] "1576737.37471207" "RAN"
## [71,] "781856.046829388" "RAN"
## [72,] "1101375.98446798" "RAN"
## [73,] "1817092.60446859" "RAN"
## [74,] "2105.24224905522" "LOWCARBON"
## [75,] "1826534.74527454" "RAN"
## [76,] "1155480.88628367" "RAN"
## [77,] "1471672.36197739" "RAN"
## [78,] "633901.745067233" "RAN"
## [79,] "911610.216606713" "RAN"
## [80,] "2290213.70950657" "RAN"
## [81,] "2646262.40850179" "RAN"
## [82,] "1352243.52082662" "RAN"
## [83,] "1265953.96448741" "RAN"
## [84,] "2100616.21162727" "RAN"
## [85,] "1776678.07569972" "RAN"
## [86,] "836.397054014272" "LOWCARBON"
## [87,] NA "CRASHED"
## [88,] "438550.74712463" "RAN"
## [89,] "1279.62367587908" "LOWCARBON"
## [90,] "1524159.38243703" "RAN"
## [91,] "2188559.27733915" "RAN"
## [92,] "1050366.8900468" "RAN"
## [93,] "1429623.48469536" "RAN"
## [94,] "27456.120078662" "LOWCARBON"
## [95,] "1398.90443438924" "LOWCARBON"
## [96,] "1037573.57180142" "RAN"
## [97,] "80867.6892480567" "LOWCARBON"
## [98,] "1744445.9170165" "RAN"
## [99,] "1233881.05138578" "RAN"
## [100,] "1443212.18749931" "RAN"
## [101,] "1898806.28032103" "RAN"
## [102,] "22.6523656647353" "LOWCARBON"
## [103,] "1618940.96335165" "RAN"
## [104,] "1773612.19063594" "RAN"
## [105,] "1505628.99582318" "RAN"
## [106,] "1511165.57724497" "RAN"
## [107,] "2387524.55902676" "RAN"
## [108,] "1808140.75254821" "RAN"
## [109,] "3034612.74359176" "RAN"
## [110,] "1305467.18461631" "RAN"
## [111,] "2215.47450906594" "LOWCARBON"
## [112,] "339510.74735257" "RAN"
## [113,] "1073457.42123026" "RAN"
## [114,] "1004856.00646112" "RAN"
## [115,] "1580083.43485615" "RAN"
## [116,] "921086.751260579" "RAN"
## [117,] "897774.501909049" "RAN"
## [118,] "2678468.32019488" "RAN"
## [119,] "975974.054951735" "RAN"
## [120,] "1304973.40986444" "RAN"
## [121,] "2795729.48708828" "RAN"
## [122,] "972523.559775463" "RAN"
## [123,] "1049919.96813186" "RAN"
## [124,] "2997.38035283134" "LOWCARBON"
## [125,] "561.565402536467" "LOWCARBON"
## [126,] "2988253.78284404" "RAN"
## [127,] "2081895.14453078" "RAN"
## [128,] "3070797.46549112" "RAN"
## [129,] "1403506.15570796" "RAN"
## [130,] "827048.517335001" "RAN"
## [131,] "395143.022994698" "RAN"
## [132,] "2163311.88964093" "RAN"
## [133,] "2178050.99534827" "RAN"
## [134,] "566076.622930143" "RAN"
## [135,] "1392961.37882784" "RAN"
## [136,] "18392.8206574553" "LOWCARBON"
## [137,] "1947440.97630319" "RAN"
## [138,] "2636818.64751266" "RAN"
## [139,] "28499.059770153" "LOWCARBON"
## [140,] "1164715.18341187" "RAN"
## [141,] "0" "LOWCARBON"
## [142,] "835144.438570216" "RAN"
## [143,] "2128420.55576015" "RAN"
## [144,] "1521308.9730314" "RAN"
## [145,] NA "CRASHED"
## [146,] "2748024.4267935" "RAN"
## [147,] "973221.738622719" "RAN"
## [148,] "860049.23757151" "RAN"
## [149,] "1324456.75170062" "RAN"
## [150,] "1201010.03107696" "RAN"
## [151,] "1142603.5225206" "RAN"
## [152,] NA "CRASHED"
## [153,] "2693611.00987939" "RAN"
## [154,] "1069541.52591967" "RAN"
## [155,] "2393741.91555386" "RAN"
## [156,] "1382026.59899902" "RAN"
## [157,] "1862826.09944603" "RAN"
## [158,] "819708.358546789" "RAN"
## [159,] "16.1599645777015" "LOWCARBON"
## [160,] "10.6491420068659" "LOWCARBON"
## [161,] "1237522.03995008" "RAN"
## [162,] "851733.274253283" "RAN"
## [163,] "1499290.2664036" "RAN"
## [164,] "969437.761077347" "RAN"
## [165,] "654971.149328721" "RAN"
## [166,] "1182209.44900622" "RAN"
## [167,] "2118636.567432" "RAN"
## [168,] "1990818.87407025" "RAN"
## [169,] "2257432.14856896" "RAN"
## [170,] "1550902.99283446" "RAN"
## [171,] NA "CRASHED"
## [172,] "8616.48535791419" "LOWCARBON"
## [173,] "2063335.56382181" "RAN"
## [174,] "2217.33780499236" "LOWCARBON"
## [175,] "1455273.15033427" "RAN"
## [176,] "1602432.8716883" "RAN"
## [177,] "1442763.37340397" "RAN"
## [178,] "935002.418487836" "RAN"
## [179,] NA "CRASHED"
## [180,] "2398227.96374005" "RAN"
## [181,] "2257472.96108099" "RAN"
## [182,] "1463089.42164426" "RAN"
## [183,] "968805.302928887" "RAN"
## [184,] "724344.993211166" "RAN"
## [185,] "202447.955948166" "RAN"
## [186,] "1436879.80366549" "RAN"
## [187,] "1050598.86322165" "RAN"
## [188,] "858303.300545437" "RAN"
## [189,] "837.8470560364" "LOWCARBON"
## [190,] "77.7376830697535" "LOWCARBON"
## [191,] "2666235.23760768" "RAN"
## [192,] "1442330.64357137" "RAN"
## [193,] "300301.613025114" "RAN"
## [194,] "1867448.54134229" "RAN"
## [195,] "2318123.19343285" "RAN"
## [196,] "1192558.07059905" "RAN"
## [197,] "296580.206341884" "RAN"
## [198,] "198948.53460339" "RAN"
## [199,] "1475476.98479445" "RAN"
## [200,] "1151974.92734534" "RAN"
## [201,] "1045780.41572662" "RAN"
## [202,] "597.349982204987" "LOWCARBON"
## [203,] "2088874.69446935" "RAN"
## [204,] "1001228.25695345" "RAN"
## [205,] "1163274.13917084" "RAN"
## [206,] NA "CRASHED"
## [207,] "16.9065379434552" "LOWCARBON"
## [208,] "1719257.89125913" "RAN"
## [209,] "1223443.31154094" "RAN"
## [210,] "436285.944826903" "RAN"
## [211,] "2517584.52687686" "RAN"
## [212,] "1147738.92352667" "RAN"
## [213,] "16630.6275934419" "LOWCARBON"
## [214,] "1136523.58134658" "RAN"
## [215,] "2945929.85082532" "RAN"
## [216,] "1027358.29145306" "RAN"
## [217,] "1090954.7293875" "RAN"
## [218,] NA "CRASHED"
## [219,] "474.410621824389" "LOWCARBON"
## [220,] "1182142.21931005" "RAN"
## [221,] "1244732.72703289" "RAN"
## [222,] "817223.408383973" "RAN"
## [223,] "789494.766979139" "RAN"
## [224,] "1286529.489797" "RAN"
## [225,] "771006.894030684" "RAN"
## [226,] "778690.42112912" "RAN"
## [227,] "177.882608699528" "LOWCARBON"
## [228,] "2084474.06576161" "RAN"
## [229,] "1225.95556037021" "LOWCARBON"
## [230,] NA "CRASHED"
## [231,] "1756439.8256892" "RAN"
## [232,] "14878.4358912422" "LOWCARBON"
## [233,] "309973.949628558" "RAN"
## [234,] "8764.0532200681" "LOWCARBON"
## [235,] "267339.183244262" "RAN"
## [236,] "5418.55514931159" "LOWCARBON"
## [237,] "1441975.92917763" "RAN"
## [238,] "709437.349215471" "RAN"
## [239,] NA "CRASHED"
## [240,] "259956.693826904" "RAN"
## [241,] "88222.8675543907" "LOWCARBON"
## [242,] "4766.6621059493" "LOWCARBON"
## [243,] "143.9677977123" "LOWCARBON"
## [244,] "1137573.83629277" "RAN"
## [245,] "655601.191706854" "RAN"
## [246,] "594181.861682584" "RAN"
## [247,] "488056.40357227" "RAN"
## [248,] "1128079.61589957" "RAN"
## [249,] "2897145.33655448" "RAN"
## [250,] "2726652.51399176" "RAN"
## [251,] "21571.315698287" "LOWCARBON"
## [252,] NA "CRASHED"
## [253,] "1961574.09889377" "RAN"
## [254,] "47946.6186543254" "LOWCARBON"
## [255,] "27076.4135515304" "LOWCARBON"
## [256,] "2251606.48396996" "RAN"
## [257,] "3273070.39797162" "RAN"
## [258,] "645886.579584998" "RAN"
## [259,] "586088.226578541" "RAN"
## [260,] "1156781.11188133" "RAN"
## [261,] NA "CRASHED"
## [262,] "1031063.75077349" "RAN"
## [263,] NA "CRASHED"
## [264,] "1917263.39494636" "RAN"
## [265,] "1615558.13924909" "RAN"
## [266,] "2533649.88926883" "RAN"
## [267,] "1531972.92447036" "RAN"
## [268,] "1509.83741632229" "LOWCARBON"
## [269,] "2231990.01459277" "RAN"
## [270,] "39017.5736788715" "LOWCARBON"
## [271,] "0.0789806248899001" "LOWCARBON"
## [272,] "2324156.32963896" "RAN"
## [273,] "21.1635773171758" "LOWCARBON"
## [274,] "588731.278489821" "RAN"
## [275,] "8434.87646631176" "LOWCARBON"
## [276,] "1052314.17562924" "RAN"
## [277,] "233938.318100738" "RAN"
## [278,] "1353008.97324927" "RAN"
## [279,] "7.79517493656395" "LOWCARBON"
## [280,] "2224267.59887657" "RAN"
## [281,] "2480640.44093978" "RAN"
## [282,] "8334.02080434849" "LOWCARBON"
## [283,] "0" "LOWCARBON"
## [284,] "2112123.26551258" "RAN"
## [285,] "3575717.57030972" "RAN"
## [286,] "1137829.11665267" "RAN"
## [287,] "768796.012526191" "RAN"
## [288,] "182756.792189772" "RAN"
## [289,] "2807794.29937303" "RAN"
## [290,] "26186.6324336018" "LOWCARBON"
## [291,] "5.16383244612062" "LOWCARBON"
## [292,] "5272823.4048439" "RAN"
## [293,] "3378889.69905921" "RAN"
## [294,] "967239.187151969" "RAN"
## [295,] "1513055.5460534" "RAN"
## [296,] "853512.111943073" "RAN"
## [297,] "1195395.23093059" "RAN"
## [298,] "1191767.22530273" "RAN"
## [299,] "2981.94763769056" "LOWCARBON"
## [300,] "1564848.16964212" "RAN"
## [301,] "2177454.63091984" "RAN"
## [302,] "1585858.82741674" "RAN"
## [303,] "1402252.58495003" "RAN"
## [304,] "1875165.36531898" "RAN"
## [305,] "1428993.72501249" "RAN"
## [306,] "2041031.91696762" "RAN"
## [307,] "1812481.75362187" "RAN"
## [308,] "1499581.42629162" "RAN"
## [309,] "1398649.12948884" "RAN"
## [310,] "1622463.15225339" "RAN"
## [311,] "2021446.30744064" "RAN"
## [312,] "1528131.98863308" "RAN"
## [313,] "1346977.76387442" "RAN"
## [314,] "1803078.20336496" "RAN"
## [315,] NA "CRASHED"
## [316,] "3600.55093489759" "LOWCARBON"
## [317,] "2735351.14292133" "RAN"
## [318,] "1875802.01944773" "RAN"
## [319,] "1011538.95542792" "RAN"
## [320,] "1571171.02786292" "RAN"
## [321,] "1441178.46552292" "RAN"
## [322,] "1449150.98896031" "RAN"
## [323,] "634697.85860513" "RAN"
## [324,] "1559508.88789955" "RAN"
## [325,] "2079199.4549005" "RAN"
## [326,] "31732.7657074611" "LOWCARBON"
## [327,] "44021.6825703995" "LOWCARBON"
## [328,] "1107846.07986755" "RAN"
## [329,] "1622331.38157197" "RAN"
## [330,] "2030849.88858994" "RAN"
## [331,] "468054.403308293" "RAN"
## [332,] "1203175.42133573" "RAN"
## [333,] "53313.7875880245" "LOWCARBON"
## [334,] "1091233.61100635" "RAN"
## [335,] "1731524.89117561" "RAN"
## [336,] NA "CRASHED"
## [337,] "2281894.39736456" "RAN"
## [338,] "1856521.22642198" "RAN"
## [339,] "21.6520858369119" "LOWCARBON"
## [340,] NA "CRASHED"
## [341,] "2285885.08156568" "RAN"
## [342,] "694644.015393941" "RAN"
## [343,] "2571.46567734245" "LOWCARBON"
## [344,] NA "CRASHED"
## [345,] "1864233.90728213" "RAN"
## [346,] "1759935.59114204" "RAN"
## [347,] "1119324.48885291" "RAN"
## [348,] "2252758.98892193" "RAN"
## [349,] "1625761.18856225" "RAN"
## [350,] "781589.125235808" "RAN"
## [351,] "1244575.03675289" "RAN"
## [352,] "2355091.83872167" "RAN"
## [353,] "675511.888141798" "RAN"
## [354,] "1098030.76399517" "RAN"
## [355,] "1273025.56057886" "RAN"
## [356,] "836188.644024318" "RAN"
## [357,] "2597925.33982223" "RAN"
## [358,] NA "CRASHED"
## [359,] "679816.768584808" "RAN"
## [360,] "1211960.97253689" "RAN"
## [361,] "143589.482061787" "RAN"
## [362,] "1087122.24802111" "RAN"
## [363,] "2382389.99621423" "RAN"
## [364,] "1809180.07632715" "RAN"
## [365,] "2249280.13648273" "RAN"
## [366,] "1882529.41130647" "RAN"
## [367,] "5328.34619164727" "LOWCARBON"
## [368,] "1208986.32211504" "RAN"
## [369,] "1107490.30053689" "RAN"
## [370,] "2216991.27203092" "RAN"
## [371,] "893859.502962013" "RAN"
## [372,] "76.2542326522061" "LOWCARBON"
## [373,] "1312674.28606182" "RAN"
## [374,] "415447.916753839" "RAN"
## [375,] "850881.513824418" "RAN"
## [376,] "885339.13594881" "RAN"
## [377,] "2113588.20086514" "RAN"
## [378,] "1425799.60275817" "RAN"
## [379,] "1128105.10960357" "RAN"
## [380,] "2371649.3296751" "RAN"
## [381,] "255886.390305109" "RAN"
## [382,] "2466945.24454644" "RAN"
## [383,] "2033112.04297303" "RAN"
## [384,] "1190784.9478732" "RAN"
## [385,] "206780.13645465" "RAN"
## [386,] "431400.751720456" "RAN"
## [387,] "1808897.43954395" "RAN"
## [388,] "221742.90645669" "RAN"
## [389,] "1499767.91545002" "RAN"
## [390,] "825569.031312678" "RAN"
## [391,] "2219524.63248659" "RAN"
## [392,] "1751328.0315877" "RAN"
## [393,] "141464.11415107" "RAN"
## [394,] "1905921.73027521" "RAN"
## [395,] NA "CRASHED"
## [396,] "86129.1602837143" "LOWCARBON"
## [397,] "1738959.2987119" "RAN"
## [398,] "2400853.04029598" "RAN"
## [399,] "595.41911312446" "LOWCARBON"
## [400,] NA "CRASHED"
## [401,] "2431040.53251091" "RAN"
## [402,] "1616335.42542376" "RAN"
## [403,] "1412281.70818034" "RAN"
## [404,] "1930697.97345133" "RAN"
## [405,] "2288.77902682165" "LOWCARBON"
## [406,] NA "CRASHED"
## [407,] NA "CRASHED"
## [408,] "5679.83400627306" "LOWCARBON"
## [409,] "2255483.7662744" "RAN"
## [410,] "1697200.48172835" "RAN"
## [411,] NA "CRASHED"
## [412,] "86772.1408945445" "LOWCARBON"
## [413,] "58.8614698072582" "LOWCARBON"
## [414,] "8079.09621326771" "LOWCARBON"
## [415,] "1222547.82839755" "RAN"
## [416,] "1643981.26262162" "RAN"
## [417,] "1090488.74884965" "RAN"
## [418,] "1844070.58227165" "RAN"
## [419,] "1435918.61989897" "RAN"
## [420,] "2331558.94159436" "RAN"
## [421,] "352790.572217638" "RAN"
## [422,] "2612299.38776824" "RAN"
## [423,] "2289592.04066901" "RAN"
## [424,] "966660.438227921" "RAN"
## [425,] "608747.910345104" "RAN"
## [426,] "3814774.23854009" "RAN"
## [427,] "25651.4463067504" "LOWCARBON"
## [428,] "859941.950259722" "RAN"
## [429,] "1175141.7966137" "RAN"
## [430,] "1388753.39016613" "RAN"
## [431,] "1307268.4862548" "RAN"
## [432,] "1230418.86711069" "RAN"
## [433,] "27562.79145666" "LOWCARBON"
## [434,] "664304.835670389" "RAN"
## [435,] "2309620.06944286" "RAN"
## [436,] "1442045.89806" "RAN"
## [437,] "1080196.27863663" "RAN"
## [438,] "12971.8130503748" "LOWCARBON"
## [439,] "2112311.00279278" "RAN"
## [440,] "1865902.63590031" "RAN"
## [441,] "1131873.97181111" "RAN"
## [442,] "2225694.00032936" "RAN"
## [443,] "2487387.00847225" "RAN"
## [444,] "18754.0635224065" "LOWCARBON"
## [445,] "153238.051869403" "RAN"
## [446,] "1900215.30615228" "RAN"
## [447,] "1413233.49798306" "RAN"
## [448,] "115.430806874736" "LOWCARBON"
## [449,] "166814.283585256" "RAN"
## [450,] "1707727.0507292" "RAN"
## [451,] "904634.041106401" "RAN"
## [452,] "126458.029982815" "RAN"
## [453,] "1192391.15413939" "RAN"
## [454,] "1167369.26995129" "RAN"
## [455,] "512992.72511025" "RAN"
## [456,] NA "CRASHED"
## [457,] "66995.4914433489" "LOWCARBON"
## [458,] "1500316.82656525" "RAN"
## [459,] "778062.924384724" "RAN"
## [460,] "19313.4679437398" "LOWCARBON"
## [461,] "1358187.39438262" "RAN"
## [462,] "2141129.84897055" "RAN"
## [463,] "7961.46211772983" "LOWCARBON"
## [464,] "1761228.27665708" "RAN"
## [465,] "667103.947334787" "RAN"
## [466,] "1516391.6244355" "RAN"
## [467,] "2069132.40213798" "RAN"
## [468,] "2188.4705746542" "LOWCARBON"
## [469,] "1156158.26689581" "RAN"
## [470,] "3456195.90521345" "RAN"
## [471,] "2280381.4764162" "RAN"
## [472,] "679642.282835393" "RAN"
## [473,] "1.07229600832304" "LOWCARBON"
## [474,] "2476588.47000487" "RAN"
## [475,] "838463.746365504" "RAN"
## [476,] NA "CRASHED"
## [477,] "1458452.40164747" "RAN"
## [478,] "44044.2655378839" "LOWCARBON"
## [479,] "651655.234920384" "RAN"
## [480,] "1970699.20064016" "RAN"
## [481,] "1405.80514674083" "LOWCARBON"
## [482,] "78.8471668056193" "LOWCARBON"
## [483,] "14637.1482404126" "LOWCARBON"
## [484,] "732601.372513896" "RAN"
## [485,] "446047.596040638" "RAN"
## [486,] "1838908.81647584" "RAN"
## [487,] "236713.773033354" "RAN"
## [488,] "2479198.89697359" "RAN"
## [489,] "1279945.11633116" "RAN"
## [490,] "56645.3158685308" "LOWCARBON"
## [491,] "3255383.23008106" "RAN"
## [492,] "597093.74530109" "RAN"
## [493,] "1529677.02897881" "RAN"
## [494,] "787714.938366425" "RAN"
## [495,] "2798.98766749579" "LOWCARBON"
## [496,] "319.650226473177" "LOWCARBON"
## [497,] "832705.134925607" "RAN"
## [498,] "2919543.18256786" "RAN"
## [499,] "513561.801092756" "RAN"